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Stimulated small-angle light scattering was used to measure the structure factor of a Bose-Einstein 
condensate in the phonon regime. The excitation strength for phonons was found to be significantly 
reduced from that of free particles, revealing the presence of correlated pair excitations and quantum 
depletion in the condensate. The Bragg resonance line strength and line shift agreed with predictions 
of a local density approximation. 
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Spectroscopic studies have been used to assemble a 
complete understanding of the structure of atoms and 
simple molecules. Similarly, neutron and light scattering 
have long been used to probe the microscopic excitations 
of liquid helium [|]-|| , and can be regarded as the spec- 
troscopy of a many-body quantum system. With the ex- 
perimental realization of gaseous Bose-Einstein conden- 
sates, the spectroscopy of this new quantum fluid has 
begun. 

The character of excitations in a weakly-interacting 
Bose-Einstein condensed gas depends on the relation be- 
tween the wavevector of the excitation q and the inverse 
healing length ^~^ = ^/2racs/h which is the wavevector 
related to the speed of Bogoliubov sound Cs = y/Ji/rn 
where ^ = Aiih ana/m is the chemical potential, a the 
scattering length, tiq the condensate density, and m the 
atomic mass. For large wavevectors {q ^ C~^)j the ex- 
citations are particlelike with a quadratic dispersion re- 
lation. Excitations in the free-particle regime have been 
accessed by near- resonant light scattering [p|. For small 
wavevectors {q <^ C~^)j the gas responds collectively 
and density perturbations propagate as phonons at the 
speed of Bogoliubov sound. Such quasi-particlc excita- 
tions have been observed at wavelengths comparable to 
the size of the trapped gas |6| and thus were strongly 
influenced by boundary conditions. 

In this Letter, we describe the use of Bragg spec- 
troscopy to probe excitations in the phonon regime. Us- 
ing two laser beams intersecting at a small angle, exci- 
tations in a Bose-Einstein condensate were created with 
wavevector q < f ~^, thus optically "imprinting" phonons 
into the gas. The momentum imparted to the condensate 
in the form of quasi-particles was directly measured by a 
time-of- flight analysis. This study is the first to explore 
phonons with wavelengths much smaller than the size of 
the trapped sample, allowing a direct connection to the 
theory of the homogeneous Bose gas. By direct compari- 
son, we show the excitation of phonons to be significantly 
weaker than that of free particles. This provides dra- 
matic evidence for the presence of correlated momentum 
excitations in the many-body condensate wavefunction. 

In optical Bragg spectroscopy, an atomic sample is il- 



lumined by two laser beams with wavevectors ki and 
k2 and a frequency difference to which is much smaller 
than their overall detuning A from an atomic resonance. 
The intersecting beams create a periodic, traveling in- 
tensity modulation /mod(r,t) = /cos(q • r — ust) where 
q — ki — k2. Atoms exposed to this intensity modu- 
lation experience a potential due to the ac Stark effect 
of strength Vmod = hV'^/SA x Iraod/hat 0, from which 
they may scatter. Here F is the linewidth of the atomic 
resonance and Igat the saturation intensity. 

The response of a many-body system to this per- 
turbation can be evaluated using Fermi's golden 
rule. We express Vmod in second-quantized notation 
Vmod = VI2 (pt(q)e-'-* -f pt(_q)e+-*) where /5t(q) = 
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dk is the Fourier transform of the atomic density 



k ^k+q 

operator at wavevector q and a^ (aj,) is the destruction 
(creation) operator for an atom with momentum ?ik. For 
the ground state \g) with energy Eg, the excitation rate 
per particle is then 
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where excited states |/) have energy Ef, N is the num- 
ber of atoms in the system, and w/?, ~ V/2h is the two- 
photon Rabi frequency. Thus, light scattering directly 
measures the dynamical structure factor, ^(q, w), which 
is the Fourier transform of density correlations in state 
\g) |||,gl. Integrating over w gives the static structure 
factor ^(q) = (5|p(q)pt(q)|5)/7V. 

In this work, measurements were performed on both 
magnetically trapped and freely expanding Bose-Einstein 
condensates of sodium. Condensates of « 10^ atoms were 
created by laser and evaporative cooling and stored in a 
cigar-shaped magnetic trap with trapping frequencies of 
ujr = 2Tr X 150 Hz and ojz — 27t x 18 Hz in the radial and 
axial directions, respectively 0]. 

The condensate was then exposed to two laser beams 
which intersected at an angle of 14° and were aligned 
symmetrically about the radial direction, so that the 
difference wavevector q was directed axially (Fig. pi). 
Both beams were derived from a common source, and 



then passed through two acousto-optical modulators op- 
erated with the desired frequency difference cj, giving 
both beams a detuning of 1.6 GHz to the red of the 
|F = 1) -^ \F' = 0,1,2) optical transitions. Thus, at 
the optical wavelength of 589 nm, the Bragg recoil ve- 
locity was hq/m = 7.2 mm/s, giving a predicted Bragg 
resonance frequency of a;° = fiq^ /2m ~ 2tt x 1.5 kHz for 
free particles. The beams were pulsed on at an intensity 
of about 1 mW/cm^ for a duration of 400 ^s. To suppress 
superradiant Rayleigh scattering flQ], both beams were 
linearly polarized in the plane defined by the condensate 
axis and the wavevector of the light. 
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FIG. 1. Observation of momentum transfer by Bragg scat- 
tering, (a) Atoms were exposed to laser beams with wavevec- 
tors ki and k2 and frequency difference lj, imparting momen- 
tum hq along the axis of the trapped condensate. Absorp- 
tion images (b,c) after 70 ms time of flight show Bragg scat- 
tered atoms distinguished from the denser unscattered cloud 
by their axial displacement. Curves (d,e) show radially aver- 
aged (vertically in image) profiles of the optical density after 
subtraction of the thermal distribution. The Bragg scattering 
velocity is smaller than the speed of sound in the condensate 
(position indicated by circle) . The response of trapped con- 
densates (b,d) was much weaker than that of condensates after 
a 5 ms free expansion (c,e). Images are 3.3 x 3.3 mm. 

The effect of Bragg scattering on a trapped condensate 
was analyzed by switching off the magnetic trap 100 /iS 
after the end of the light pulse, and allowing the cloud 
to freely evolve for 70 ms. During the subsequent expan- 
sion, as the density of the atomic cloud dropped, quasi- 
particles in the condensate transformed into free particles 
and were then imaged by resonant absorption imaging 
perpendicular to the magnetic trap axis (Fig y). Bragg 
scattered atoms were distinguished from the unscattered 
atoms by their axial displacement. The speed of Bogoli- 
ubov sound at the center of the trapped condensate is re- 
lated to the velocity of radial expansion Vr as Cg = v^j \f2 
ijlll (ca = \\ mm/s at [ijh = 6.7 kHz as shown in Fig. |^). 
Thus, comparing the axial displacement of the scattered 
atoms to the radial extent of the expanded condensate, 
one directly sees that the Bragg scattering recoil velocity 
is smaller than the speed of sound in the trapped conden- 
sate, i.e. the excitation in the trapped condensate occurs 
in the phonon regime. 

For comparison, Bragg scattering of free particles was 
studied by allowing the atomic sample to freely expand 
for 5 ms before application of a light pulse at equal inten- 



sity [|2|. During the 5 ms expansion, the atomic density 
was reduced by a factor of 23 and the speed of sound by 
a factor of 5 from that of the trapped condensate. Thus, 
Bragg scattering in the expanded sample occurred in the 
free-particle regime. 

The momentum transferred to the atomic sample was 
determined by the average axial position in the time-of- 
flight images. To extract small momentum transfers, the 
images were first fitted (in regions where the Bragg scat- 
tered atoms were absent) to a bimodal distribution which 
correctly describes the free expansion of a condensate in 
the Thomas-Fermi regime, and of a thermal component 
||l3| . The chemical potential /i of the trapped condensate 
was determined from the radial width of the condensate 
distribution [gTj . The non-condensate distribution deter- 
mined by the fit was subtracted from the images before 
evaluating the momentum transfer. 

By varying the frequency difference lj, the Bragg scat- 
tering spectrum was obtained for trapped and for freely 
expanding condensates (Fig. ^. The momentum transfer 
per atom, shown in units of the recoil momentum hq, is 
anti-symmetric about w = as atoms are Bragg scattered 
in either the forward or backward direction, depending 
on the sign of w [n4| . 
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FIG. 2. Bragg scattering of phonons and of free particles. 
Momentum transfer per particle, in units of Tig, is shown vs. 
the frequency difference tj/27r between the two Bragg beams. 
Open symbols represent the phonon excitation spectrum for 
a trapped condensate at a chemical potential /x//i = 9.2 
kHz. Closed symbols show the free-particle response of an 
expanded cloud. Lines are fits to the difference of two Gaus- 
sian line shapes representing excitation in the forward and 
backward directions. 

From these spectra, we determined the total line 
strength and the center frequency (Fig. ||) by fitting the 
momentum transfer to the difference of two Gaussian line 
shapes, representing excitation in the forward and back- 
ward direction. Since S'(q) = 1 for free particles, we 
obtain the static structure factor as the ratio of the line 
strengths for the trapped and the expanded atomic sam- 
ples. Spectra were taken for trapped condensates at three 
different densities by compressing or decompressing the 
condensates in the magnetic trap prior to the optical ex- 
citation. 
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Substituting the Bogoliubov operators, one obtains ||T^ 
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FIG. 3. (a) Static structure factor S(q) and (b) shift of 
the hue center from the free particle resonance, ^(q) is the 
ratio of the line strength at a given chemical potential n to 
that observed for free particles. As fi increases, the structure 
factor is reduced, and the Bragg resonance is shifted upward 
in frequency. Solid lines are predictions of a local-density ap- 
proximation using ujq = 27rx 1.38 kHz. Dotted line indicates 
a mean-field shift of 4/i/7/i as measured in the free-particle 
regime (data from M shown in open symbols). 

The Bragg resonance for the expanded cloud was cen- 
tered at 1.54(15) kHz, in agreement with the expected 1.5 
kHz response for free particles [ [l5[ , with an rms- width of 
900 Hz consistent with Doppler broadening |1^ . The re- 
sponse of trapped condensates was strikingly different. 
As the density of the trapped condensates was increased, 
the Bragg scattering resonance was significantly weaker 
in strength and shifted upwards in frequency. This re- 
flects the changing character of the excitations created 
by Bragg scattering as the speed of sound was increased: 
at a fixed Bragg scattering momentum, the excitations 
passed from the free-particle to the phonon regime. 

To account for this behavior, we use the Bogoliubov 
description of the interacting homogeneous Bose-Einstein 
condensate [0. The hamiltonian 
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for a gas in volume V, where huj^ = h k'^/2m, is ap- 
proximated by replacing the zero-momentum operators 
with c-numbers aj = Sq = \/iVo where A'o is the num- 
ber of atoms with zero momentum. Neglecting terms 
of order N^^^^, the hamiltonian is diagonalized by a 
canonical transformation to operators defined by dk — 
Ukbk — Vkh__j, where Uk = coshc/)^, Vk = sinh^fc and 
tanh20fe — n/{hLU^. + /i). The energy of the Bogoliubov 
excitation created by S[, is hcjj^ = yj huj^ifioj^ + 2jjl). 

Neglecting small contributions representing multiparti- 
cle excitations ||,^ , the single quasi-particle contribution 
to the static structure factor is 

'5'(q) = -^{g\{aqa\ + d^_^a^q + a^_^d\ + dqd^^\g). (2) 
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5(q) ~ (U^ + Vf^ - 2UqVq) = U^ U ^ 
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In the limit fnu^ 3> /x, the Bogoliubov excitations become 



identical to free-particle excitations {uq 
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and S'(q) ^ 1. For phonons (7iw° <C /i), S'(q) — > 
hq/2mcs ~ q^/\/2 and the line strength diminishes lin- 
early with q. 

To the same order of approximation, the quasi-particle 
resonance is undamped, and the dynamical structure fac- 
tor is simply S'(q, cj) = S{q)S{uj — oj^) (satisfying the 
/-sum rule: J ujS{q^,uj)duj — cj^ ||^). Thus, accompany- 
ing the diminished line strength, the Bragg resonance 
is shifted upward from the free particle resonance by 

Equivalently, the suppression of the Bragg resonance 
in the phonon regime can be understood in terms of 
the many-body condensate wavefunction. The static 
structure factor is the magnitude of the state vector 
1^) ~ '^k^k+n^k\g) I \N . The macroscopic population 
of the zero-momentum state picks out two relevant terms 
in the summation: 
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These represent two means by which momentum is im- 
parted to the condensate: either by promoting a zero- 
momentum particle to momentum ?iq, or else by demot- 
ing a particle from momentum — ?iq to zero momentum. 

If correlations could be neglected, the total rate of ex- 
citation would simply be the sum of the independent 
rates for these two processes, proportional to (e"'"|e+) = 
{Nq) + 1 = w^ and {e-\e-) = (N-q) = v^ where (Nk) is 
the expected number of atoms of momentum Kk in the 
condensate. This would apply, for example, to a conden- 
sate in a pure number state, or to an ideal gas condensate 
with a thermal admixture of atoms with momenta ±?iq, 
and would always lead to 5(q) > 1. 

Yet, for the many-body ground state of the interacting 
Bose gas, the behavior is dramatically different. Colli- 
sions of zero-momentum atoms admix into the conden- 
sate pairs of atoms at momenta ±?iq, the population of 
which comprises the quantum depletion [|l9| . As a re- 
sult, the two momentum transfer mechanisms described 
above produce indistinguishable states, and the rate of 
momentum transfer is given by the interference of two 
amplitudes, not by the sum of two rates. Pair excita- 
tions in the condensate are correlated so as to minimize 
the total energy, and thereby give destructive interfer- 
ence between the two momentum transfer processes, i.e. 
S'(q) — [uq — VqY < 1. For high momentum, {Nq) <^ 1 
and the interference plays a minor role. In the phonon 
regime, while the independent rates ui and Vq (and hence 
{N±q)) diverge as 1/q, the highly correlated quantum de- 
pletion extinguishes the rate of Bragg excitation. 



These results for the homogeneous Bose gas can be di- 
rectly applied to trapped, inhomogeneous condensates by 
a local density approximation since the reduced phonon 
wavelength q~^ (0.4 /im) is much smaller than the size 
of the condensate (> 20 /im) [p0 21|. In the Thomas- 
Fermi regime, the condensate has a normalized density 
distribution f{n) — 15n/4no-\/l — n/riQ where ng is the 
maximum condensate density. The Bragg excitation line 
shape is then 
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which can be integrated to obtain the line strength 5(q) 
and center frequency. The line strength has the limit- 
ing values of S'(q) — > 157r/32 (?ia;°/2/i)^/^ in the phonon 
regime and ^(q) ^ 1 — Aji/Thuj^ in the free-particle 
regime [^21. In accordance with the /-sum rule, the cen- 
ter frequency uj is given as a;^/5(q). 

These predictions are shown in Fig. ^ using loq — 2'kx 
1.38 kHz which accounts for the expected mean-field shift 
of the Bragg resonance for the expanded condensate. 
Both the line strength and the shift of the Bragg res- 
onance are well described by our treatment. For compar- 
ison, previous measurements ||5|] of the mean-field shift of 
the Bragg resonance (Afi/Th) in the free-particle regime 
are also shown, clearly indicating the many-body charac- 
ter of low energy excitations. 

In conclusion, stimulated light scattering was used 
to excite phonons in trapped Bose-Einstein condensates 
with wavelengths much smaller than the size of the 
trapped sample. The static structure factor was shown 
to be substantially reduced in the phonon regime. This 
modification of light-atom interactions arises from the 
presence of a correlated admixture of momentum exci- 
tations in the condensate. The observed reduction of 
S'(q) also implies a reduction of inelastic Rayleigh scat- 
tering of light with wavevector fc by a condensate when 
hu>^ < ^ [^. This effect may reduce heating in opti- 
cal dipole traps and reduce the optical density probed in 
absorption imaging. For example, the absorption of near- 
resonant light by a homogeneous sodium condensate at 
a density of 3 x 10^^ cm^'^ |24| should be reduced by a 
factor of two. 
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